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Abstract 

Imagine that you have several sets of two coupled qubits, but you do not know the parameters of 
their Hamitonians. How to determine these without resorting to the usual spectroscopy approach 
to the problem? Based on numerical modeling, we show that all the parameters of a system of 
two coupled qubits can be determined by applying to it an external classical noise and analysing 
the Fourier spectrum of the elements of the system's density matrix. In particular, the interlevel 
spacings as well as the strength and sign of qubit-qubit coupling can be determined this way. 
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INTRODUCTION 



Despite steady successes in fabrication and measurement techniques, the experimental 
characterization of multi-qubit systems jl, 2] remains a challenge due to their complicated 
level structure. Our goal here is to determine the system's parameters, as distinct from the 
more dimcultproblem of determining its state, which has to be tackled using quantum state 
tomography 3|,|4|. For example, neither the strength nor the sign of the qubit-qubit coupling 
are known a priori. One of several standard approaches studies the resonant response of 
quantum macroscopic systems to an external coherent signal (see, e.g., |5-7|), allowing to 
determine the qubit parameters by scanning the frequency range of the external signal. The 
difficulty in the straightforward application of this approach, due to the fact that only few 
qubits can be actually accessed, and the relation of this problem to the general field of 
inverse problems, were addressed in 

An alternative approach to the standard spectroscopic methods of characterization would 
use as a drive a broad-band noise. We call it active noise spectroscopy, as distinct from the 
"passive" noise spectroscopy of Ref. [5], where the response of the noise spectrum to a 
coherent monochromatic drive was measured. 

Recently, we have shown [lOj that classical noise applied to a qubit produces persistent os- 
cillations of the off-diagonal density matrix elements ( "coherences" ) despite finite dephasing 
and relaxation times. In other words, a moderate amount of external noise enhances quan- 
tum coherence, which manifests in oscillations with a frequency corresponding to quantum 
transitions between the ground and first excited states. There exists an optimal noise am- 
plitude: at lower noise level, oscillations are suppressed, while as the noise is increased, the 
oscillations become random and the corresponding spectroscopic peak is eventually smeared 
away. Indeed, for zero noise, the oscillations of the off-diagonal elements of the density 
matrix decay on the time scale of r, where 1/r is the dephasing rate. Moderate phase- 
insensitive noise excites the system from time to time, allowing the qubit to evolve with its 
own frequency between the relatively rare noise spikes, thus, uncovering quantum dynam- 
ics. Strong noise produces strong spikes very often, thus leaving no time for the coherent 
evolution. This phenomenon is related to both classical and quantum stochastic resonances, 
which manifest in various physical systems (see, e.g., llMl6|). 

In this paper we investigate how these effects of classical noise can help determine the 
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parameters of a multiqubit system. Specifically, we consider two coupled qubits and analyze 
the spectrum of the density matrix excited by white Gaussian classical noise. We numerically 
show that the resulting noise spectra contain four peaks, which correspond to the interlevel 
transitions in the system. From these, the energy spectrum and all the model parameters 
of the qubits are readily obtained. In addition, the correlations in the matrix elements 
corresponding to different qubits can be used to conclude whether the qubits are coupled 
ferro- or antiferromagnetically. 



MODEL 



Two coupled qubits can be described by the Hamiltonian [17] 



H=-\Y. [V* • <Ancr> • gotf (1) 

where a J z and a 3 x are Pauli matrices corresponding to either the first (j = 1) or the second (j = 
2) qubits, and the eigenstates of a{ are the basis states in the localized representation of the 
jth qubit at zero coupling. Note that the results obtained below do not qualitatively depend 
on the type of coupling (e.g., cr],a1 versus cr^crjj): in any case the noise will allow to determine 
the parameters of the two-qubit Hamiltonian. For this reason, and for demonstrating the 
physical principles of noise-induced spectroscopy, we consider two identical qubits. The 
tunneling splitting energies Ai^ (in case of the identical qubits we will be investigating here: 
Ai = A2 = A) are determined by the design and fabrication details of the device, while the 
bias energies 6j(t) can be controlled externally and, in our case, are only driven by the noise, 

e J (t) = ^(t). (2) 

The Gaussian white noise considered here is zero-averaged and delta-correlated: 

($(*)) = 0, (UmAt)) = 1D5 hr 5(t - 0- (3) 

where D is the noise intensity, which should be defined for each particular system (see, e.g., 
the example of two flux qubits described below). The uncorrelated noise sources affecting 
the qubits ( "local" noise) tend to be more detrimental to their quantum coherence than the 



correlated ones 



22], which makes Eq. (j3J the "worst case scenario". 
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Master equation 



By writing the qubit density matrix p as 

. 1 



a,b=0,x,y,z 



we can rewrite the master equation 



in the form 



n 02/ 
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~dl 
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H(t),p + fp 
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-ei(t)n W 3 - e 2 (t)n ZJ/ - (ri + r 2 )(n 22 - z T1 z T2 ) 



Here we used the standard approximation for the dissipation operator f via the dephasing 
and relaxation rates to characterize the intrinsic noise in the system. Also, hereafter we 



assume for simplicity that relaxation rates are the same for both identical qubits, i.e., T<pi = 
r«£2 = and T r i = r r2 = r r , and that the temperature is low enough, resulting in 
the equilibrium values of the diagonal elements of the qubit density matrices being Z? 2 = 
Z T1 = 1. All the simplifying assumptions (e.g., Ai = A 2 , T rl = T r2 , = etc.) do not 
qualitatively affect our results reported below. For instance, if A x ^ A 2 , the spectrum in 
Fig. 2 will have more peaks, corresponding to larger numbers of levels due to the lifting of 
the artificial degeneracy. 

In the limit of zero coupling (g = 0), there exists a solution of Eqs. ([5]) with no en- 
tanglement between qubits. This solution can be written as a direct product of two 
single-qubit density matrices written through the corresponding Bloch vectors: pj = 
|(1 + Xjf x + YjT y + Zjf z ). The components of what can be called the Bloch tensor U a b 
are then all zero except for (Jl ox , Il oy , Il oz ) = (X u Y 1 , Z x ) and (IL^, IT^, Ii zo ) = (X 2 , Y 2 , Z 2 ). 
If the interaction is not zero, the entanglement between these qubits generates all the com- 



ponents of the Bloch tensor to be non-zero 22 



23] and such an entangled state persists on 



the time scale 1/r after the interaction is later switched off [g(t > to) = 0]. 

This reflects the fact that, in the presence of interactions, the eigenstates of the system are 
. and the noise terms in the eigenbasis will thus maintain the off-diagonal 



entangled 



22 



terms in the density matrix of the two-qubit system. 



Two flux qubits 

As a specific example of our approach, which can be experimentally implemented, we pro- 
pose to measure two (almost) identical superconducting flux qubits consisting of a supercon- 
ducting loop interrupted by four Josephson junctions and coupled via a coupler loop [24j (See 
Figure 1). The state of each qubit is controlled by the applied magnetic flux $w) = / e ^$ 
through the loop, where $o is the flux quantum. In the vicinity of f^ 1 ' = fjp = 1/2, the 
ground state of the system is a symmetric superposition of the states \L) and \R), with a 
clock- and counterclockwise circulating superconducting current I p , respectively. In the basis 
{\L}, \R)} the two-qubit system can be described by the Hamiltonian (TjQ) with ej = I p &o5f^ 
with classical flux bias fluctuations 5f e in the qubit loops around 1/2, while the tunneling 
amplitude A is determined by the fabrication of the loop and the junctions. Note that the 
components of the density matrix can be measured directly, e.g., by monitoring the current 
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fluctuations in the flux qubits: Ij(t) = I p X(t). The direct relation of this spectrum to the 
current /volt age noise spectrum in the resonant tank (LC) circuit coupled to the qubit was 
used in Ref. 



SIMULATION RESULTS 

Using the dimensionless time i = tA, we numerically solved the system (jSJ) by the Ito 
method for two coupled qubits driven only by white classical noise, choosing parameters for 
damping T^/A = r r /A = 0.1 close to the ones experimentally found in flux qubits. The 
spectra of X\ = H ox and Z\ = H oz , for g = 0.5, are shown in Fig. 2. Since this two-qubit 
system is only driven by noise, the spectrum of both X\ and Z\ is enhanced by increasing the 
noise, and peaks become more distinguished if noise is not too high. These spectra exhibit 
four maxima, whose positions nicely agree with the frequencies of the interlevel transitions 
(in units of A): 



1-kvx = cut = 2g, 2nis 2i 3 = w 2 , 3 = yl + g 2 ± g, 2-kv^ = u; 4 = 2y 1 + g 2 (6) 

which have values v\ ~ 0.16, v 2 ~ 0.1, z/3 rj 0.26, and v± ~ 0.36. Two peaks out of these 
four frequencies are clearly seen on the Sx spectra in Fig. 2, while the other two peaks are 
better seen on the Sz spectra. Either of these two spectra is sufficient to measure both the 
coupling constant g and the tunneling splitting energy A, while the remaining spectrum can 



be used for control. Note here that, unlike the single-qubit case [10|, there are peaks on 
both Sx and Sz even for small values of the coupling strength g, which illustrates our earlier 
remark on the entangled nature of the eigenstates revealed by classical noise. 

To determine whether the coupling is "ferro-" or "antiferromagnetic" , that is, the sign 
of the coupling constant g, we study the time correlations in the density matrix elements 
tl x(t) — Xi{t) an d n xo (t) = X 2 (t), for g = ±0.7 and g = 0. Numerically solving equations 
(J5]) we obtained the time sequences Xj(ti) shown in Fig. 3, where ti is the discretized time 
of the simulation. Correlations and anticorrelations are clearly seen for ferromagnetic and 
antiferromagnetic coupled qubits, while almost no correlations are seen for the decoupled 
ones. A qualitative physical picture of these correlations in the time domain is readily 
understood. For instance, for ferromagnetic coupling, the Bloch vectors of the two qubits 
tend to allign for weak enough noise. A stronger noise excites partially-coherent oscillations 
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in the intervals between two sequential noise spikes, but the qubit-qubit oscillations still 
tend to preserve the ferromagnetic ordering (which results in the correlations seen in Fig. 
3) even for the dynamicaly evolving qubits. Similarly the antiferromagnetic coupling tends 
to produce antt-correlations in the qubit dynamics, as seen in Fig. 3. 

To quantitatively describe these correlations we plot the sample Pearson correlation co- 
efficient 

r = I 2 I 2 

y/nEilPM - {^UU)) ^n^KAU) - (EMU)) 2 
as a function of the coupling constant g (bottom panel of Fig. 3). Here n is the total number 

of simulation time steps. The module of the correlation coefficient r exhibits a maximum at 

\g\ 0.7. At larger \g\ the oscillations become weaker, since the uncorrelated external noises 

in two qubits suppress each other via their coupling, so that the noise-induced oscillations 

weaken. The sign of r coincides with the sign of the coupling g, which allows to easily 

distinguish between ferro- and antiferromagnetic couplings. 

CONCLUSIONS 

We have demonstrated that quantum correlations in a two-qubit system can be high- 
lighted by the presence of classical noise. As an application of this effect, we suggest the use 
of noise spectroscopy. Namely, the measurement of the fluctuation spectra of the system, as 
a means to determine the relevant parameters of the multiqubit system. 
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FIG. 2: (Color online.) Spectral density Sx{w) (top panel) and Sz{oj) for two values of the noise 
(D/A = 0.04 and D/A = 0.013) and normalized coupling g/A = 0.5. Four peaks, two per panel, 
can be easily distinguished, and these correspond to the four inter level frequencies v\, 1/2,^3, and 
1/4. The insets show their corresponding time sequences X(t) and Z(t). 
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FIG. 3: (Color online.) Time sequences for X\{t) = Tl ox {t) (continuous red curve) and X%{t) = 
n xo (i) (dot-dashed black curve) for values of the coupling constant g = ±0.7 (top two panels) and 
(third panel). The anticorrelations (top panel, g > 0) and correlations (second panel, g < 0) 
are clearly seen for nonzero coupling (|<7| = 0.7), while there are no correlations for g = 0. The 
bottom panel shows the dependence of the correlation coefficient r on the coupling constant g, 
when D/A = 0.013. 



11 



